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(54) IMPROVEMENTS IN OR RELATING TO SEISMIC PROSPECTING 

/-7i\ w,. rwFVRON RESEARCH COMPANY, a corporation duly organised under 
the lals of ih^ S^aYe of DdS-^^^^^ States of America, of 525 Market Street. San 

Frar^c^co California United States of America, do hereby declare the mvent.on. for which 
we JS thaTa patent mav be granted to us. and the method by which .t .s to be performed. ^ 
S tr, hp mrticularlv described in and bv the followmg statement.- 

ThisTnS on>eSs to the art of seismic prospecting for petroleum reservoirs by 
mnUiole"DTnt survevine techniques, and more particularly to the art of converting 
E^oh PnteK reflS^on^ amph anomalies associated with one or more common 
Serpoin^J obSSd on seismic record traces into diagnostic indicators of the presence of 
10 gas in the underlying =^"^bs"^f/;^J;^^ ^or petroleum has involved the creation of 
acrust,rdiYturbances"bov^^^^^^^^ the surface of the earth, using explosives, 

air^ur^s or laree mechanical vU^rators' Resulting acoustic waves propagate downwardly m 
orth .nH nr^nT.rtiiliv rctlected back toward the surface when acoustic impedance 

paSlaV bver depends on the velocity and density content between that layer and the 
layer which overlies it. say according to the formula 
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C Reflect = AR/Ai = Vzd. - Vid, 

V,d^ + V.d, 
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25 where AR is the amplitude from the reflected sigrial and Ai is the amplitude of the incident 
<;ian il- V, is thc velocitv of the wave in the overlying medium 1. V2 is the \ciocit> in inc 
rlildtL laver beloS contact line; d, is the density of thc ovcrlymg medium 1 : and d, is 

'^\^S^^r:'^^^^c':::il.c.c^ acoustic waves were recorded immediately in 3^ 
the field as visible stde-bv-side. dark, wigglv lines on white paper (' seismograms ). At 
p'res'enl the iS reproductions - in a digital format - °" "-f^'^^^^P^^ ^IJI^.^iSg 
are reduced to visible side-by-side traces on paper or film in large central computing 

Ar'sti'ch centers sophisticated processing makes possible thc distinction of signals from 
no^e in cases that Sild have seemed homeless in the early days of seismic pr"spectmg_ 
Until 1965 almost all seismic surveys conducted used an automatic gam control which 
continuous V adjusted the aa.n of amplifiers in the field to account for decreasing amounts 
oreSv rom late reflection arrivals. As a result, reflection coefficients could not be 
accu^atllv determined However, with the advent of the expander circuit and binary gam 
amSuf e s et n of the amplifiers can now be controlled and amplitudes recorded precisely; 
Ts makes Pt possible to conserve not only the special characteristics of thc reflections, but 

'''^^iS'S^lu^'p^^en. in the prior equipment. Cornputers of^^^^^^^^^^^ 
thc use of a comparison technique because ot their small word size and tin> core storage 
45 Today more powerful computers with array processors and economical floating point 
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capabilities now enable modern geophysicists to maintain control of the amplitude of all 
recorded signals. The '*floating point" capability is especially effective in expanding 
computer work size by a large factor and in eliminating the need for computer automatic 
gain control. 

5 " That is to say, in summary, as a result of the above advances, reflections from many 5 
thousands of feet below the earth's surface can now be confidently detected and followed 
through sometimes hundreds of side-by-side traces, the shortening or lengthening of their 
corresponding times of arrival being indicative of the shallowing" or deepening of actual 
sedimentary strata of interest. Still, as a general rule, all that can be hoped for the seismic 

10 reflection method is to detect stratigraphic interfaces and the interfaces as they deviate from 10 
horizontality of these interfaces, so that subsurface "structures" could be defined in which 
oil or gas might possibly be trapped. 

Apropos of the above has been use of ultra-high amplitude anomalies in seismic traces to 
infer the presence of natural gas in situ. Seismic interpreters have used so-called 

15 '^bright-spot " analysis to indicate several large gas reservoirs in the world, especially in the 15 
Gulf Coast of the United States. Such analysis is now rather common in the oil industry, but 
it is not without its critics. Not only cannot the persistence of such increased amplitude 
anomalies be taken as confirmation of the lateral extent of the gas reservoir, but also the 
anomaly itself (in some cases) may not represent reflections of a discontinuity of a 

20 gas-bearing medium and its over- or underlying associated rock strata. E.g.. experience has 20 
shown that in certain situations, similar phenomena occur which can confuse the 
interpreter. E.g.. if the shape of the horizon is such that it focuses the energy back to the 
surface, it may increase the amplitude of one or more of the records akin to reflections from 
gas-saturated strata. Lithology of the horizon - singly and in combination - can also have a 

25 similar effect, producing high-amplitude reflections "in the absence of gas within the pore 25 
space of the stratum of interest. Examples of the latter: conglomeratic zones, hard streaks 
of silt or lime and lignite beds. 

The present invention improves the ability of the seismologist to correctly differentiate 
high-intensity anomalies of multiplc-point-coverage seismic traces of gas-bearing strata 

30 from those of similarly patterned reflections of other types of stratigraphic configurations 30 
containing no gas accumulations. 

An object of the invention is the provision of a novel method of correctly differentiating 
high-intensity anomalies provided by multiple-point seismic traces of gas-bearing structures 
from those of a similarly patterned intensity associated with strata contaimng no gas 

35 accumulation. 35 
Fn accordance with the present invention, interpretation of high-intensity seismic events 
from traces obtained from multiple-point coverage of a subterranean earth formation using 
an array of source means and detectors adjacent to the earth's surface is obtained to 
indicate gas-bearing strata in a highly surprising and accurate manner. After the field data 

40 have been obtained in which the data of common centcrpoints are associated with more 40 
than one source-detector pair, the data arc indexed (*^iddressed'') whereby all recorded 
traces are indicated as being a product of respective source-detector pairs of known 
horizontal offset and centerpoint location. Thereafter, high-intensity amplitude anomalies 
in said traces are correctly associated with gas-bearing strata on a surprisinglv accurate 

45 selection basis: amplitude intensity of said anomaly rnust change -- progressively in an 45 
increasing or decreasing manner - as a function of horizontal offset. A furlhcr refinement 
of the method of the present invention may be in order under some circumstances. A single 
common-centcrpoini trace, or even single common-ccntcrpoint gathers, may have a major 
drawback in such cases -- poor signal-to-noisc ratios. As a result, progressive amplitude 

50 change as a function of offset cannot be resolved. In accordance with this'invention and as a 50 
means of signal enhancement, trace summations can prove beneficial in improving record 
resolution, say on a basis of a stacking "window" having a two-dimensional index for 
addressing the traces: X common offset^dimension long by Y common centerpoints wide. 
For example, where 2400^;r> common-ccnterpoint stacked traces have been obtained (i.e., 
55 24 traces per gather) by multiple-point-coveragc field techniques, each gather can in turn be 55 
"de-stacked" to provide original but corrected locaiional traces. Then on the ba.sis of a 
stacking window four (4) common offset dimensions long bv five (5) common centerpoints 
wide, several such traces, say 10, can be stacked and^the stacked trace displayed as a 
function of offset. 

60 Result: progressive change in amplitude intensity as a function of similar intensity changes 60 
in offset can be more easily observed whereby in situ gas is, more likelv than not. in the 
pore space of the structure of interest. 
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^'^^ plan S o?a ^Td of cemcrpoints produced in ^^-^^^.^.^^ 5 

P^^^,"^"!'Tmoaefortvp?^^^ reflecting horizons within an earth formation that can be 
associated wuh the charac^^^ ,0 
seiSic w ve?l-;ciSe'd'wit°h honzons of Fj.ure 2 wh.ch a.d .n the 

germination ot .the prese^^^^^^^^^^^^ mathematical nature. 

asSction if pe ccm4e^^^^^^^ iHustratin'g the relationship of Po.sson s rat.o 

^^HilSrcpi^rofS^^^^^^ 

wh^'fem" geo"ltdc"l'?ia'n\t^rmaToThas oc'curred to better illustrate processes assoc.atcd 

'''%;^s9^^<^(i^lr!n^^^^^ akin to those shown in F.gures 7 and 

20 8 frcarrv^ng out the method of the present invention, usmg a programmed digual 20 

'°F/g»r^fo andTi' are schematic diagrams of elements within the digital computing system 

"^Ss n-vire true seismic record sections and portions of sections, illustrating the 
dia Josdc cl.p"abmtv"rthe method of the present invention in pred.ctmg the presence of 25 

^^^Beir discussion ofan emhSment of the invention within an actual field environment^ 
a brief description of the mathematical and theoretical concepts behind it ma> prove 

''pr/s^^v ' if mav be Kr'St to'?ndicate lithologv limitations associated with the present 30 

Poisson-s ratio and resulting high-intensity amplitude anomalies provided on seismic faces. 
35 While Poisson's ratio (o) has the general formula 
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45 where V„ is compressional velocity and \\ is shear velocity of the med.un . the ^^^cpt iocs 
have Dhvsical sicnificancc. For example, consider a slender cvhndr.cal rod ot an elastic 
matcrSSappI a compressional force to the ends. .As the rod cMianges shape, the length 
of he rod will decrease bv /IL. while the radius will increase byAR. Poisson s ratio is 
^fSed L the ado ot the r-elative change in radius (^R/R) to the relative change in Icng h 

50 Hence a compres.sible material has a low Poisson s ratio, while an incompressible 50 

TTTitpri'il a liauid) has a hitih Poisson s ratio. 

S equation above also indicates the relationship of the compressional and shear wave 
veiocitiel of the material. V„ and V. respectively, re., that .'^'^'fff " ^o;";'"^ ' ^. 

determined dvnamicallv bv measuring the P-wave and S-wave velocities. Onlv two ot th. 

^'^IL^nrt^hrhed'^uEs' on- r'c^^i^ctTo^n and transmission seismic waves useful in 

'7r)''Koifo%'d''o".°r95.s"tr'"On the Effect of Poisson s Ratios of Rock Strata in the 
Reflection Coefficients of Plane Waves". Geophysical P'-^^P^^ting^ V«L ... r^^^^ 4. 
60 r) Koefoed O.. 1962. for 'Reflection and Transmission Coetticients tor l lane 
\ nnairiidinal Incident Waves ". Geophvsical Prospecting. Vol. 10. No. .v 

AT Muskat N^^ Meres. M.W.*. 1940. for -Reflection and Transmission Coett.cients 
for Pl'ine Waves in Elastic Media". Geophvsics. Vol. No. ^. 

(4> Toolcv R.D. Spencer. T.W. and Sagoa H.F.. for • Reflection and Transmission of 
65 Plane Compressional Waves". Geophysics. Vol. 30. No. 4 (196?). 
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(5) Cosiain, J.K., Cook, K.L. and Algcrmisshi, S.T., for "Amplitude, Enersv and 
Phase Angles of Plane SP Waves and Their Application to Earth Crustal Studies'V Bull. 
Seis. Soc. Am.. 53. p. 1639 et. scg. 

All of the above have focused on'^the complex modeling of variation in reflection and 
transmission coefficients as a function of angle of incidence. 5 

The problem is complicated, however. E.g".. isoiropfc media with laver index of the 
strata. i=l for the incident medium and i=2 for the underlying medium, have been 
modeled using equations for P-wave reflection coefficient A^, and for P-vvave transmission 
displacement amplitude coefficient Apj. In such equations, the value of Poisson's ratio (o) is 
required, since both P-wavc and S-wave velocities arc utilized. For each of the media, i.e.. 10 
the incident or underlying medium, three independent variables exist: P-\vave velocity, o 
and bulk density, or a total of six variables for both media. But for a sincle interface, only 
four independent variables were required: p-wave velocity ratio, the densuv ratio. Poisson's 
ratio in the incident medium and Poisson s ratio in the underlying medium. 
15 But to provide for the many combinations of possible variations.'the above-listed studies 15 
have either: 

(a) generated many (literally thousands) plots of a mathematical nature for various 
parameters, values in which there was little relationship with true gcophvsical applications, 
since the latter were hopelessly obscured and unappreciated; or 
20 (b) made simplistic assumptions that, although using actual calculations, nevertheless 20 
did not express the true nature of transmission and reflection coefficients, in particular 
hthological situations associated with the accumulation of gaseous hydrocarbons within an 
actual earth formation. 

In summary, while reference (2) concluded that chanee in Poisson's ratio in the two 
25 bounding media can cause change in the reflection coefficient as a function of an^le of 25 
incidence. (2) did not relate that occurrence to litholoev associated with the accumulation 
of gaseous hydrocarbons in the surprising manner of the present invention. 

The present invention teaches that gas-containing strata have low Poisson's ratios and 
that the contrast with the overburden rock as a function of horizontal offset produces a 
30 surprising result; such contrast provides for a significant - and prouressive - chance in 30 
P-wave reflection coefficient at the interface of interest as a function of^angle of incidence of 
the incident wave. Thus, differentiating between high-intensitv amplitude anomalies of 
nongas and gas media is simplified by relating progressive chanize'in amplitude intensity as a 
function of offset between source-detcctor^pairs, i.e.. ansle^of incidence beine directly 
35 related to offset. " " 35 

Also, the behavior of P-wave travel as a function of litholosv and horizontal offset 
between a respective source-detector pair associated with a t;ivcn"iocational trace provide 
the following amplitude response signatures of interest: 

d) where the gas-containing media are gas sands underlvinii shale, such as found in the 
40 Gulf Coast, amplitude responses increase with offset: ^ 40 
(2) where the gas-containing media include limestone undcrlvint: shale, such as found in 
the North Sea, the amplitude anomalies of the interface decrease with offset. 

Now m more detail, attention should be directed to the Fiizures. particularly Fiuure 1. 
Note that, inter alia. Figure 1 illustrates in some detail how ithe terms of intc'rcst^in this 
45 application arc derived: e.g.. the term "centcrpoint'' is a tzeoeraphical location located 45 
midway between a series of sources S,, S.-.-S^ of a geophysical field svsicm 9 and a set of 
detectors D,. D2. .D„, at a datum horizon near the earth's surface, the centerpoints are 
designated C,, Q.-.Cp in the Figure, and are associated with a trace derived bv placement 
of a source at that centerpoini location followed immediately by relocatiniz a detector 
50 thereat. ' ' " 

I.e.. if the sources S,...Sn are excited in sequence at the source locations indicated, traces 
received at the different detector locations shown can be related to common centerpoints 
therebetween. If such traces arc summed, a gather or iiroup of traces is formed. I.e., if the 
reflecting interface is a flat horizon, the depth point where reflection occurs will define a 

55 vertical line which passes through the ccnterpoint of interest. Applvinu static and dynamic 
corrections to the field traces is equivalent (under the above facts) to pTacinc the individual 
sources S.-.-S,, at the ccnterpoint in sequence followed bv replacement with the 
detectors D,...Dn, ot interest at the same locations. If the traces associated with a common 
ccnterpoint are summed, a series of enhanced traces, sometimes called CDPS (Common 

60 Depth Point Stack) traces, is provided. 

Figure 2 illustrates reflection phenomena of a three-lavcr mode! typical of a youn<^ 
shallow geologic section 10. such as found in the Gulf Coast, illustrating how reflecti(m 
phenomena associated with the traces associated with the field system 9 of Figure I can he 
related to the presence of gas. 

65 Section 10 includes a gas slmd 1 1 embedded in a shale stratum 12. Assume a Poisson's 65 
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ratio of 0 1 for the eas sand and of 0.4 for the shale, a 20% velocity reduction at interface 
13 say from 10 OOO'/sec to SOOOVsec. and a 10% density reduction from 2.40 g/cc to 2.16 

The actual P-wave reflection coefficient Apr can be related to section 10 by Equation (1 ) 
S be7ow- also^ P wave tra^^^^^ displacement amplitude coefficient Apt can s.m.larh he 5 
related in accordance with Equation (2) below. 



10 A =L^^^-^ (1) 
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h,-2aiki(cLV + c,i) 
'^P' h,(./' + T + X) 



^ = kt k? 6(a,C2 + a^c,) (3) 



10 



(2) 

15 



20 


= b-yy + a,c,v- 

X 








20 


25 


= a,c,(5' ~ 4a,c, |- b^) 

r\ = 6e2 — El 


(5) 
(6) 






25 


30 


g = 6 - I 


(7) 






in 


35 


V = Se, 2b- 
1 = E, + 2f.b- 
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b = h|Sin6 
at = h{ - b^ 
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Q = kt - b^ 
hi = 1/Vpi 
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pi = 


S.wave velocity 
Defisity 
Layer Index 


55 


60 


ki = 1/V,i 


(17) 


e = 


Angle of Incidence 
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Equations (1) and (2) are, of course, the two basic equations of wave travel in an earth 
formation and are for isotropic media with the layer index being i = l for the incident 
medium and i = 2 for the underlying medium. Equations (3) through (17) simpiv define 
intermediate variables. 

5 As an example of calculations associated therewith, if 9=0" (normal incidence), the 5 
P-wave reflection coefficient Apr is equal to about -0.16 and +0.16. respectivelv. 

Figure 3 illustrates change in reflection coefficient as a function of angle of incidence 0 for 
the three-layer model of Figure 2. 

Note that solid lines 20, 21 illustrate the effects of reflection (and transmission, by 
10 omission) on the top and base of the gas sand. In line 20, at 9=0°, note that the Ap^ equals 10 
-0.16; while at 9=40^ the A^, is about -0.2S. That is, rather a surprisingly large change in 
the reflection coefficient as a function of angle of incidence occurs, with the greatest chanee 
occurring between 6 = 20° and 6=40"*. 

For the bottom layer, line 21 changes at about the same rate, but in opposite sisn. I.e. . at 
15 9=(P. Ap, is about +(J.16 and at 9=40^ Ap^ is about +0.26. Again, the sreatest\^hanse in 15 
Ap, occurs between B=20^ and B=40\ As a result, the amplitude oflhe seismic wave 
reflected from this model would increase about 70% over the angle of incidence range 
shown, i.e.. over the incremental 40 degrees shown. " ° 

While angles of incidence equal to 40^"'may seem a little large for reflection profiling 
20 (heretofore, most data arriving beyond 30° being thought useless and muted out)", 20 
experience has now ncvertheless^shown that reflection data can and do arrive at reflection 
angles greater than 30^. Hence, the angles of incidence must be determined, and one of the 
more important techniques, the staight-rav approach to estimate such angles of incidence 
(using depth-to-reflector and shot-to-detector and-shot-to-iiroup offset)^ is as set forth 
25 below: " 25 

Bi = arc tan (X/2Z) (18). 

where 6, is the angle of incidence; X is the shoi-io-dctector or shot-to-group offset and Z is 
30 the retlector depth. Velocity changes with depth can likewise be" accommodated by 30 
assuming section velocity change is o"f the form V,=V;,+KZ where K is a constant so that 
all ray paths arc arcs of circles having centers V„/K above the reference plane of interest, 
say the earth's surface. Thus, the approach should be in accordance with 

35 Bi = arc tan ZX + V.,X/K 35 

Z- + 2VJZ/K)-(X-'^) (19) 



Having now established a firm mathematical and theoretical basis for the process of the 
present invention, perhaps a description of how a geological section containinu no gas ^0 
therein would affect impedance contrast is in orderT Figure 4 illustrates the chanses^in 
reflection coefficient as a function of angle of incidence Bin the manner of Fi^ure 3. but in 
which the gas sand 10 of Figure 2 contains no gas, simulating, e.g.. a low-velocity, 
brine-sauirated. young sandstone embedded in shale. 
"♦5 The solid lines 22. 23. representing reflection coefficients, are seen to be about horizontal ^5 
between 0=0'' and 40'. slightly decreasing in magnitude as the angle of incidence increases, 
i.e., as B approaches 40°. In the above example.^il should be noted that the Poisson's ratio 
of the sandstone was assumed to he 0.4. 

Figure 5 illustrates yet another plot associated with a three-layer model akin to that 
shown m Figure 2. but in which the sandstone contains iias but is' buried deep below the 50 
earth's surface. The values for the three-laver model of Fiuurc 2 are asain used except that 
the velocity change from shale to sand is only 10%. or from lO.OOOVscc to 9(H)0'/scc. As 
shown, curves 25, 26 are even more significant: both curves are seen to increa.sc in 
magnitude trom over the 40° of change in the anizle of incidence. However, field results 
have not ventied these results, since Poisson's ratio in such gas sands may be strongly 55 
aftccted by depth, and not be as low as is now surprisingly taught bv the present invention! 

Figures 6(a), 6(b). 6(c) and 6(d) offer a possible expfanatiTm for low Poisson's ratio in 
gas-contaming strata in general and in gas sands in particular. In the Ficures. various 
quantities arc plotted as a function of percentaee of saturation. In Fiaure 6(a), P-wave 
velocity IS so plotted: m Figure 6(b), S-wave velocity Is depicted: in Figure 6(c). the ratio of 60 
Vp/Vs IS the value of interest: and in Figure 6(d), Poisson^s ratio is shown as a function of 
percent gas saturation. 

Note that Figures 6(a) and 6(b) are for sandstones buried at 6000 feet with 35% porosiiv. 
Figures 6(c) and 6(d) result from Figures 6(a) and 6(b) usins appropriate equations. Hut in 
Figure 6(d). Poisson's ratio drops from about 0.3 to 0.1 from 0% to 10% "as saturation- on 65 
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the other hand, the same ratio changes very little from 1^ to 10U9c gas saturation 

■"Henfe ^r'^^m' the abovl mathematical and theoretical concepts, displays of reflcct.on data 
Hence, iruin mcu . reflection coefficient as a function ot angle ol 

can now used to indicate ^^^^^^^^^ hydrocarbons. Such data are also now 

^ l-oni'^nfentlv avai labfe sav ESna todaVs conventional field-gathering techniques mvolvmg 
^irTn^^ areS ccfve^ ce since h? former can be derived from and is compatible with one o 
multiple-area coverage since me reflection seismic data: common-depth-point 

;CD-p)%Tth;rs"" Ani ^to^res ive chants in reflection amplitude, vs. shot-to-de.ector 
in arouD) oftse can Sm he basis of such a determination, since otlsel of any parucula 
lf;u;ce-deteaor pair iT directly related to the angle of incidence m accordance with 

^'^Bnt' mrreaea iSonal trace data, before stacking in accordance with CDPS techniques 
often hTve piuM^n^^ ratios. Thus, changes in amplitude vs. otlset may be d.ll.cul. 

'%tn/cl ?s"a'dtLri:m which illustrates a data ' addressing- technique which improves 
..mnRudc versu o^^^^ resolution in such situations: in the Figure the traces were 
SSeS usSn end-shooting array of 48 detectors with source and detectors advancing 
one detector interval per shot point -n 

on Result- a 74-fold CDP-stacked record section was generated. 

Sotffu ther each centerpoint is associated with 24 separate traces ot .varying offset 
In order to Sometrical v associate each generated locational trace with its common 
centerpoint address auidai^cc. as provided by Figure 7. is important. To Jrstand the 
natire^of Fiuure 7. assume that the sources S , .S,. . .S„ are sequentially 'o^^« '^'l^ "J.^'P"!^'^ 25 

25 SP SP sK at the top of the Figure. Assume also that the detectors are placed in hnc w,th 25 
^h,' .o.f;ces 1 e alon" the same line of survev A at the detector locations D,.p2...D,n. 
After c4h source ifact vated. reflections are received at the detectors, at the locations 

between SP, and each of the detectors, at D,. D^.-D is a series of '-'^"'f^'^ P.*f 
DyuNccn ^ , associated with a trace. In this regard and tor a turlhcr 

data" can be related to the common centerpoints midway between individual source pomis 
nnd detectors as discussed in the above-noted reference. 
40 But bv sS a field technique, data provided generate 24 separate traces assicated w,th 40 
the same centemc!^nt C,...Q. In «rde[ to index ( -address ") these traces as a function o 
sSei^iKactor^lncludinu horizontal offset and centerpoint location, a stacking chart 44 as 

^'"^har;"44';f: dia. ';am 'in whS'f t^ace is located along a plurality of oblique comnrion 
45 profile ini PL, PC:.... between a series of common offset and centerpoint locations a y( 45 
de-rees to eLh other. For best illustration, focus on a single shoipo.nt. say SP,. and on a 
Jnol^de^cc or spread havinc detectors D,.D:...D,„ of Figure 8 along survey line A 
TsfLfa source Kcated at%hotpoint SP, and activated thereafter. The detector spread 
^nTsource are ''oiled" forward along survey line A ,n the direction B. l^-^'"^^^d,y:'"";1 
50 station per activation Then after detection has occurred, and it the resulting centerpoint 
pattern is olated 45= about angle 46 to profile line PL, and projected bdow the spread as m 
F nure S as a function of common off.set values and centerpoint positions, the chait 44 of 
F "U c 8 esults Of course, each centerpoint has an amplitude vs. time trace associated 
thirewith and for didactic purposes that trace can be said to project along a line normal to 

'"t^iruld^be^^mrhasized that the centerpoints provided in Figure.s 7 and S are 
geoWaphicallv located along the line of survey A in line with the source points SP, S^ 
As the locational traces are generated, the chart 44 aids in keeping a tag on each icsulting 
uai As he de ector spread and sources are rolled forward one station and the techn.qi c 

60 repeated. toTherl^^^ traces is generated associated with -■"'-\^POints on ne^^ pro i e 
line PL,. That is. although the centerpoints arc geographically ' ^ ''^,^o^Yf^ ic^^^^^^ 
position; alonu the survev^ line A of Figure 7. by rotation ''lo^j;, 'h'^. f. i^^.J' ^T., e - 
ccntcrooint pa'ttern C,\C,'...C„' can be horizontally and vertically a igned uith tcntci 
SSrev^iouslv "cne ated. l.e at common offset values (in horizontal alignnient) ceiiai. 

65 ^rp^ofnts Signed, viz. centerpoint C. aligns with C,' as shown: further C : is aligned 
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with C2'. etc. Also, there are traces that have common centerpoints. I.e., at common 
centerpoints (in vertical alignment) cenierpoint C2 aligns with centerpoini C,'. and 
cenierpoints €3X2' and C," are similarlv aligned. Thus, via chart 44, each trace associated 
with a centerpoint can be easily '"addressed" as to: 
5 (i) its actual geographical location (i.e. , along phantom lines normal to diagonal profile 5 
lines PLi^PLt... along common location lines LLj.LL:...). so that its actual field location is 
likewise easily known: 

(ii) its association with other traces along common horizontal offset lines COL,, COL: 
...COL^: and 

10 (iii) its association with still other traces alon? common vertical centerpoint location 10 
lines CPLj.CPL..... 

Also, ''addressing" the traces by (ii) and (iii) allows such traces to be easily combined 
(summed) by calling out '^vindows" within the chart in which any traces within the window 
can be summed. E.g., it has been found convenient to establish a standard window ''width'' 

15 equal to an increased group cenierpoint line (ACPL) value of say 5, and a window ''height" 15 
equal to an incremental common group offset line (ACOL) value of say 4: hence by 
indexing the intersecting window intervals on a sequential basis, summation of traces 
therein can occur. The results are summed traces which are outputted to a display on a 
side-by-sidc basis, say as a function of amplitude intensity as a function of increasing or 

20 decreasing offset between respective source-detector pair's. Actual offset values are^not 20 
required, since relative values are usually sufficient for most diagnostic purposes. 

In carrying out the above summation process on a highspeed basis, a fully programmed 
digital computer can be useful. But electromechanical systems well known' in the art can 
also be used. In either case, the field traces must first undergo static and dynamic correction 

25 before the traces can be displayed as a function of offset to determine their potential as a 25 
gas reservoir. Such correction techniques are well known in the art — see, e.g., U.S. Patent 
2,838,743. of O.A. Fredriksson et aL for **Normal Moveout Correction with Common 
Drive for Recording Medium and Recorder and/or Reproducing Means"', assigned to the 
assignee of the present application, in which a mechanical device and method are depicted. 

30 Modern processing today uses properly programmed digital computers for that task in 30 
which the data words arc indexed as a function of, inter alia, amplitude, time, datum 
height, geographical location, group offset, velocity, and are manipulated to correct for the 
angular and horizontal offset; Tn this latter environment, see U.S. 3.731.269, Judson et al. 
issued May 1, 1973, for "Static Corrections for Seismic Traces by Cross-Correlation 

35 Method'', a computer-implemented program of the above type also assigned to the assignee 35 
of the present invention. Electromechanical sorting and stacking equipment is also Well 
known in the art and is of the oldest ways of cancelling noise. Sec. for e.xample, the 
following patents assigned to the assignee of the present mvention which contain sorting 
and stacking techniques, including beam steering techniques: 
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offset value and the results displayed, the window ''width ' is incremented a common group 
ccnterpoint interval (ACPL) and the process repeated and the results displayedT 

Figure 10 illustrates particular elements of a computing system for carrying out the steps 
of Ffgures 9 (a) and 9(b). While many computing syste'ms are available to carry out the 
5 proce'ss of the invention, perhaps to best illustrate operations at the lowest cost per 5 
instruction, a microcomputing system 50 is didactically best and is presented in detail 
below. The system 50 of Figure 10 can be implemented on hardware provided bv many 
different manufacturers, and for this purpose, elements provided by Intel Corporation, 
Santa Clara. California, may be preferred. 

10 Such a system 50 can include a CPU 51 controlled by a control unit 52. Two memory units 10 
53 and 54 connect to the CPU 51 through BUS 55. Program memorv unit 53 stores 
instructions for directing the activities of the CPU 51 while data memory'unit 54 contains 
data (as data words) rcFated to the seismic data provided by the field acquisition system. 
Since the seismic traces contain large amounts of bit data, an auxiliary memory unit 55 can 

15 be provided. The CPU 51 can rapidly access data stored through addressing the particular 15 
input port, say at 56 in the Figure. Additional input ports can "also be provided to receive 
additional information as required from usual external equipment well known in the art, 
e.g., floppy disks, paper-tape readers, etc., including such equipment interfaced throuch 
input interface port 57 tied to a keyboard unit 58 for such devices. Usine clock inputs. 

20 control circuity 52 maintains the proper sequence of events required for any processing 20 
task. After an instruction is fetched and decoded, the control circuitrv issues the 
appropriate signals (to units both internal and external) for initiating the proper processing 
action. Often the control circuitry will be capable of responding to e'xlernal sicnals. such as 
an interrupt or wait request. An interrupt request will cause" the control circuitry 52 to 

25 temporarily interrupt main program execution, jump to a .special routine to service the 25 
interrupting device, then automatically return to the main proeram. A wait request is often 
issued by memory units 53 or 54 or an I/O element that operates slower than the CPU. 

For outputting information, the system 50 can include a printer unit 59 wherebv the 
amplitude of the summed traces as a function of time is printable. Of more use as an output 

30 unit, how^ever. is disk unit 60. which can temporarily store the data. Thereafter, an off-line 30 
digital plotter capable of generating a side-by-side display is used in conjunction with the 
data on the disk unit 60. Such plotters arc available in the art. and one proprietarv model 
that I am familiar with uses a computer-controlled CRT for optically mcriiing onto 
photographic paper, as a display mechanism, the seismic data. Briefly, in such a plotter the 

35 seismic data, after summation, are converted to CRT deflection signals; the resultinc beam 35 
is drawn on the face of the CRT and the optically merged record of the event indicated, say 
via photographic film. After a predetermined number of side-bv-side lines have been 
drawn, the film is processed in a photography laboratory and hard' copies returned to the 
interpreters for their review. 

40 Figure 11 illustrates CPU 51 and control unit 52 in more detail. 40 
As shown, the CPU 51 includes an array of registers generally indicated at 62 tied to an 
ALU 63 through an internal data bus 64 under control of control unit 52. The registers 62 
are temporary storage areas. Program counter 65 and instruction recistcr 66 have dedicated 
uses; the other registers, such as accumulator 67. have more general uses. 

45 The accumulator 67 usually stores one of the seismic operands to be manipulated bv the 45 
ALU 63. E.g., in the summation of traces, the instruction may direct the ALU 63 to not 
only add in sequence the contents of the temporary registers containinu predetermined 
trace amplitudes together with an amplitude value in the accumulator, but also store the 
result in the accumulator itself. Hence, the accumulator 67 operates as both a source 

50 (operand) and a destination (result) register. The additional registers of the array 62 are 50 
useful in manipulation of seismic data, since thev eliminate the need to shuffle results back 
and forth between the external memory units of Figure H) and accumulator 67. In practice 
most ALU's also provide other built-in' functions, includinu hardware subtraction, boolean 
logic operations, and shift capabilities. The ALU 63 also can utilize flau bits generated by 

55 FF unit 73 which specify certain conditions that arise in the course of arithmetical and 55 
logical manipulations. Flags typically include carrv. zero, sign, and parity. It is possible to 
program jumps which are conditiona'lly dependent'on the status of one or'more fkms. Thus, 
for example, the program may be designed to jump to a special routine if the carrv^bit is set 
tollowmg an addition instruction. 
60 Instructions making up the program for operations involving seismic data are stored in 60 
the program memory unit 53 of the CPU 5 1 of Figure 1 1 , The program is operated upon in a 
sequential manner except when instructions in the memorv units 53. 54 call for special 
commands such as '"jump" (or "call*') instructions. While the program associated with the 
P^^M^P^. '"^'cntion IS a relatively straightforward one. hence avoiding most *M"mp" and 
65 'call instructions, "call" instructions for subroutines are common m the processing of 65 



seismic data and could be utilized, if desired^ In :'cair •-tructio- hjs.a 
special way of handling subroutines of.^^^ '° foTit incre^^^^^^ 

is noS^h"e memory unU .nown as ^ 

' '^^CpS-s have different wa>. of ^^^t^^S^rSS^t u^a rSS S^^'o^f 
storase of return addresses built into ^^e CPU itselt. utrie^^ ^.^^^^ 
external memory as the stack and s.mplvrn^^ a P"^^"^'^,/^",^,^^,^ entry. The stack 
register 70, Figure 1 1 , which contains th^e address ot tne mos completed. 10 

10 thus saves the address of the 'ns^uction to be e^^^^'^J.^f;^^ Us program counter 65. The 
Then the CPU 51 loads the address specified in the ^^^^^^ Lbroutine. The last 

ne.xt instruction fetched will therefore be the first step^^^^ ^^^^^.^ 

'"Srg-n'JlA^^ien^dTsSe'^ ^he^SS^atio^ns^f fhVSu 51. Table 1 Ts presented below 
15 containing a full instruction set for its operations. 
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Examples 

Diagnostic capability provided by the method of the present invention is better illustrated 
in the Examples set forth below. 

5 Example I 5 
Seismic data were obtained over a gas field near Sacramento. California. These data, in 
CDP-stacked form, arc shown in Figure 12. The field. di.scovered in 1972, consists of a 
100-foot sand which is almost fully gas-saturated. The discovery w-ell is located at about 
SP-86 of Figure 12, with the currently developed portion of the field extending from about 
10 SP-75 to SP-115. Gas occurs at a depth of about 7000 feet, which corresponds to a time of 10 
about 1.75 seconds on the plot. 

Common-depth-point gathers from 3 locations. A, B and C of Figure 12. are shown in 
Figures 13(a), 13(b) and 13(c). Both single-fold and 10-fold summed gathers are shown for 
locations A and B. while only the summed gather is shown for location C. Shot-to-group 
15 offset for all gathers increases to the left, as indicated with the minimum and maximum 15 
trace offset distances annotated. These distances change on the summed gathers because 
the summing is done over 4 offsets. 

Note the strong amplitude increase with increasing offset at locations A and B. The 
10-fold summing obviously improves signal-to-noise ratios and an amplitude increase by a 
20 factor of about three is indicated from near to far offset. Gathers at location C, however, 20 
show no indication of amplitude increase with offset, and in fact show a decrease. This 
possibly indicates an absence of gas in the vicinity of location C. This possibility is also 
supported by the presence of a gas-water contact in a well structurally projected at about 
SP-120. 

25 25 

Example II 

Seismic data were obtained in the Fallon Basin of Nevada and arc depicted in 
CDP-gathered format in Figure 14. A well was drilled at SP-I27 in Figure 14. A seismic 
amplitude anomaly is indicated at location A at about 1.6 seconds. Upon drilling, the 
30 amplitude anomaly was found to originate from two basaltic layers, 100 feet and 60 feet in 30 
thickness. As its structural position indicates, this well was a stratigraphic test in an 
undrilled basin. 

The common-depth-point gathers at the well location are shown in Figure 15. Here, there 
is a strong indication of reflection amplitude decrease with increasing offset. This finding is 
35 consistent with the absence of gas in the geologic section and the expected Poisson's ratios 35 
for sediments and basalt. 

Example III 

Seismic data obtained from an area in the Sacramento Valley, California, are depicted in 
40 Figure 16. A well was drilled at SP-61. Note the amplitude anomaly extending from about 40 
SP-45 to about SP-90 at 1.5 seconds. However, the amplitude anomaly was found to 
originate from a high-velocity conglomerate layer. 

Shown in Figures 17(a) and 17(b) are the single-fold common-depth-point gathers at two 
locations: location A at the well and location B. about 1/2 mile to the west. The gathers at 
45 location A do indeed indicate the absence of gas, i.e., no noticeable increase in reflector 45 
amplitude with offset. However, the gathers at location B do show a slight increase in 
amplitude with offset, i.e.. possible gas. 

Example IV 

50 Seismic data were obtained for another area and are depicted in Figure IS. The possible 50 

gas-related amplitude anomalies are located (i) between SP-270 and -310 at about 1.3 

seconds and (ii) between SP-250 and -300 at about 1.0 second. 

The ten-fold CDP gathers at locations A and B of Figure 18 arc shown in Figures 19(a) 

and 19(b), respectively. Here, there do indeed appear to be indications of^'ampliiude 
55 increase with offset. In Figure 19a, the anomaly appears over a region where amplitude 55 

increases with offset. In Figure 19b, the anomalv at 1.0 .seconds is thought to be related to 

low-velocity shale 

Example V 

60 Seismic data were obtained for another area and are depicted in Figure 20. The geologic 60 
section was limestone embedded in shale. The gas-related anomalv is located over C\\q 
indicated rectangular area of the Figure. 

Here, note that for this lithology, gas is indicated by decreases in amplitude with offset, 
as shown in Figure 21 representing Cd? gathers at surface locations 102 and 103 of Figure 

65 20, as viewed respectiveiv from right to left in Figure 21. 65 



17 



1 594 633 



17 



10 



15 



25 



30 



10 



The method of the present invention as described provides a geophysicist )^';th a strong 

circumstances as known to those skilled m the art. 
^"J^iShod^of de!e!mminc the location of strata containing gaseous hydrocarhons 

compnsmg the steps of ^ record of sicnals from acoustic discontinuities 

1- ^ • • I r.o/-h rif QiiH corrected traces is ident nccl m its reiauonsnip lu 

(i) selectina a first series of indexed traces withm a common ottsct. common 
centernoint window of predetermined dimensions: and 

? 0 ^surSmTnTsaid fir'^t series of selected traces to form a summed trace. 
\ The method of Claim 2 with the additional substeps ot. ^. . i . 

35 (iii) ^nc;Smtnr the window at least in the common offset d,mens>on ,o select a 35 

''nv")'' ILmmiL' u"'seco^d series of selected traces to form a second summed trace. 

4^ The m^hod of Sm 3 in which step (d) is a side-by:side d'^P o said summed 
tr-Tres as a function of pro^ressivelv chansing composite horizontal oitset values whereby 
40 progressi4 chanes Tn sakt high-intensity, amplitude ovem of step (e) arc more easily 40 

'"f^Thi method of Claim 2 in which the dimension of the window of step ( i) is four offset 

^^!r" A Sethol of dc"rmrninL SVocatlro/^tr^rixontaining gaseous hydrocarbons 
45 which tncSe" ■lUtoScally converting an original multitrace [I,'^"^ ^5 

imorovcd section having increased resolution as to the nature ol high-intensity . mpinuoe 
evems Sated ?o e l^^ from subsurface strata possibly containing gaseous hydrocar- 
bons 'saldlmpfoved section being composed of a plurality o^^^^^Zl^!^:'^^^ 
offset-and-time traces, said original record consisting ot a plurality ot "^"'i''';'^,'- .7'^' 
SO traces o ;uT^plitude-versus-hori£ontal coordinatc-and-time. each ot said traces constituting 5O 
enercv denved in associate with a particular source-detector pair ot known horizonta^ 
S\md of known centerpoint location, and representing, in part event reflections from 

said subsurface strata, said method comprising ^'"^P*^' v .iucs ind common 
(n\ • classifvin" said original traces on the basis of horizontal offset values anu common 
55 cemerpoint Seasons, whe'reby each trace is identified by a centerpoint location common .0 55 

at least another trace and a known horizontal offset value; .,„_H-,icd with said 

(b) displavino ;,t least said each trace and said another trace '^ssocialeu u n saiu 

common cSerpoint location, as a function of progressively changing horizontal offset 

vnhie*; tn form at least a segment of said improved section, anu , . . 
60 c) dcter^ninVthe location of strata containing gaseous hydrocarbons from observ^- 60 

tion of pmarTssive^hange in a high-intensity amplitude ^.^TuXT'tormiSg s^d 

and said another trace as a function of progressive change m "'*"r^f,^"^^^ 

segment of said improved section, which indicates that said « 

acoustic impedances associated with strata containing S«s-^""\f^Vdrocarbo^^^ .another f.^ 
65 7 The method of Claim 6 in which said common centerpoint location ot said another 65 
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classified trace of step (a) is also common to a plurality of additional other traces, each 
having known but progressively changing horizontal values with respect to said another 
trace. 

8. The method of Claim 7 in which the step (b) is a side-by-side display of said each 

5 trace, said another trace and said additional other traces as a function of progressively 5 
changing horizontal offset value to form said improved section whereby progressive change 
in said Tiigh-intcnsity event of step (c) is more easily identifiable. 

9. The^ method of Claim 7 in which step (a) is further characterized by: 

(i) indexing all of said traces in two dimensions whereby each of said traces is identified 

10 in its relationship to neighboring traces on the basis of progressive changes in horizontal 10 
offset value versus progressive changes in common centerpoini location; 

(ii) selecting a first series of indexed traces within a common offset, common 
ccnterpoint window of predetermined absolute dimensions; and 

(iii) summing the first series of selected traces to form a first summed trace. 

15 10. The method of Claim 9 with the additional steps of: 15 

(iv) incrementing the window in at least the common offset dimension to select a second 
series of traces: and 

(v) summing the second series of selected traces to form a second summed trace. 

11. The method of Claim 9 in which step (b) is a side-by-side display of said summed 

20 traces as a function of progressively changing composite horizontal values to form said 20 
improved section whereby progressive change in said high-intensity amplitude event of step 
(c) is more easily identifiable. 

12. The method of Claim 9 in which the window of step (ii) has dimension of four offset 
values wide by five centerpoint locational points long. 

25 13. A method of .seismic prospecting, substantially as hereinbefore described with 25 
reference to the accompanying drawings. 
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